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Abstract—A novel PM spherical actuator based on three- 
dimensional pole array is proposed and developed in this paper. 
Conventionally, 2D pole arrays are widely employed in the design 
of spherical actuators, which constraints the system torque output 
greatly. The concept of 3D pole array is aimed at improving 
the torque performance. The torque has been analyzed and the 
corresponding analytical model is established based on curve 
fitting method (CFM) due to the importantance to real-time 
control. Magnetic force has been studied in a similar way. The 
results shows that the modeling method has a relatively high 
precision and can be further used in the real-time control. 


I. INTRODUCTION 


Devlopment of robotics, automated manufacturing process- 
es and military and aeronautics has led to strong demand 
for the devices that can achieve multiple degrees of freedom 
motions in a single joint. These devices denote as spherical 
actuators, which have drawn plenty of researchers’ attention 
for several decades due to the compact structure, high preci- 
sion, fast response and no backlash and singularities compared 
with the devices with multiple single-axis motors [1], [2]. 
The first spherical actuator was invented by Williams and 
Laithwaite in 1955 [3]. Lee et al. proposed the concept of 
spherical actuator based on the variable-reluctance principle 
and conducted researches on its magnetic field and torque [4], 
[5] in 1980s. After that, more in-depth designs and analysis [6], 
[7] were conducted. Chirikjian et al. developed a PM spherical 
stepper actuator with 12-inch diameter [8]. The corresponding 
kinematic design and the commutation problem are analyzed. 
Wang et al. conducted researches on magnetic field analysis 
and torque calculation based on the prototype of the PM 
spherical stepper motor implemented by Chirikjian et al. [9], 
[10]. Motion control algorithm and kinematic characteristics 
are analyzed furthered. Wang et al. developed permanent 
magnet (PM) spherical actuators which can achieve either 2- 
DOF motion or 3-DOF motion [11], [12]. Serious of systematic 
work has been down on the magnetic field, torque calculation 
and control schemes. Studies on spherical actuators can be also 
found in [13]-[21]. 


Torque and force characteristics are extremely important 
to the analysis of the whole actuator. Thus, the torque and 
force performances are focused in this paper. A lot of relative 


work has been conducted on the calculation of the torque of 
the spherical actuators based on the prototypes. The torque 
model in [22] and [10] are formulated based on Maxwell 
Stress Tensor (MST) method. The disadvantage of MST is 
that the result is very sensitive to discretization density and 
integration contour position [15]. The coenergy method has 
been utilized by Lee et al. [7] to derive the torque model 
of variable-reluctance spherical motor (VRSM) successfully. 
However, this method becomes unsuitable for some cases with 
large air-gap or air-core coils. Alternatively, Wang [11] and 
Yan [23] calculated the torque using Lorentz force law, which 
is suitable for stators with air-core coils and rotors with PM 
poles and can be used for real-time control. However, the 
control accuracy may not be that high due to the approximation 
during the modeling process. 


To improve the torque performance, a novel spherical 
actuator based on 3D pole arrays has been proposed in our 
previous work [24], [25]. Magnetic field has been analyzed 
in analytical, numerical and experimental ways. In this paper, 
a curve fitting method (CFM) will be adopted to build the 
torque and magnetic force model to improve the modeling 
accuracy and facilitate real-time control. This method becomes 
especially useful for the control of a special spherical actua- 
tor.And the research of magnetic force model in this paper 
will be of significance to the further design and improvement. 
The rest of this paper is organized as follows. In Section 
II, the concept and working principle of the novel spherical 
actuator are introduced briefly. In Section II, the torque model 
is established. And in section IV, the magnetic force model is 
discussed. 


II. CONCEPT DESIGN AND OPERATING PRINCIPLE 


The concept of 3D pole array patterns can be illustrated in 
Fig. 1. In traditional design of spherical actuators, the coils and 
PM poles are always mounted on 2D spherical surfaces. Those 
pole array patterns can be defined as 2D pole array patterns. 
While in the novel array patterns design, the poles of the rotor 
and stator are distributed in 3D spaces. In this case, the inside 
volume of the actuator can be utilized more effectively through 
reasonable poles arrangement. The coils with current passing 
through interact with the enhanced magnetic field produced 
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by the 3D PM poles. As a result, higher force/torque can be 
produced than the the 2D pole array spherical actuator. Fig.1 
shows one typical 3D pole pattern design. As the sketch map 
illustrates, two spherical-shelf-like rotors mounted with PM 
pole arrays are fixed together through the output shaft and two 
stators mounted with coils are connected through a shaft fixed 
on the base. The two stators and two rotors gear to each other 
alternatively. To facilitate the observation, the outer shelves of 
the spherical stator and rotor are removed from Fig. 1(a) and 
(b). Two layers of poles in radial direction are employed in 
this design. 


Pole Holder 


Coil Poles 


PM Poles 


(a) Rotor poles distribution (b) stator poles distribution 


Fig. 1. Concept of three-dimensional pole arrays: (a) Rotor poles distribution 
(b) Stator poles distribution 


To facilitate prototype manufacturing and analysis, a 3D 
pole array spherical actuator with two layers of PM poles 
and one layer of coils in radial direction is developed in the 
following study. The working principle of this machine is 
illustrated in Fig. 2. The rare-earth PM poles with alternate 
polarities are mounted along the rotor equator to produce high 
flux density around the stator. The air-core coils are assembled 
on the stator symmetric about the stator equatorial plane. By 
energizing pairs of coils in two longitudinal directions, the 
rotor can tilt in two orthogonal directions as shown in Fig. 
2 (a) and (b). Energizing all circumferential coils, the rotor 
can spin about its own axis Fig. 2 (c). Thus, by regulating the 
value and direction of input currents of in the coils, desirable 
3-DOF rotary motion can be achieved within the workspace. 
Fig. 2 (d) is the assembled actuator CAD model. 


III. FORMULATION OF OUTPUT TORQUE 


Torque model is extremely important for the real-time 
control of the spherical actuator. By utilizing the torque model, 
the required currents can be calculated from the desired 
actuator torque. Analytical torque model based on Lorentz 
force law can reflect relationship the between the torque 
performance and the actuator parameters very well and be 
applied in optimization design. However, it is not that suitable 
for the control of the spherical actuator because that too 
many approximation of the complex boundary conditions will 
deteriorate the control accuracy. The finite-element simulation 
can make up for the deficiency and be more accurate to obtain 
the torque data although it can not reflect the function relation 
between the torque performance and the actuator parameters. 
Thus, to build a more accurate torque model, it can be analyzed 
using the FEM. The FEM model is established in Maxwell 3D. 
Herein the torque is modeled through curve fit of the data from 
FEM simulation [6], [26]. 


Outer rotor 


(b) Second tilting motion 


Outer rotor 


(c) Spinning motion 


(d) Assembled actuator 


Fig. 2. 3-DOF motion of the spherical actuator. (a) First tilting motion. (b) 
Second tilting motion. (c) Spinning motion. (d) Assembled actuator 


A. Torque Analysis with FEM 


As the coils of the spherical actuator are mounted on non- 
ferromagnetic cores, the torque output is proportional to the 
current input. Therefore, the torque model can be derived from 
one coil pole model through linear superposition. The air gap 
between the inner rotor and the stator is 0.5mm. And the air 
gap between the outer rotor and the stator is 2mm. Considering 
the need of fixing the coils, the latter is a little bigger than the 
former one. In the future research, more means of coil fix could 
be discussed to diminish the air-gap’s influence on the torque 
performance. The key parameters of the rotor and stator are 
given in Table I. 


TABLE I. PARAMETERS OF THE INNER ROTOR 
Inner Rotor Cylindrical PM Value 
Outer Radius of Outer Rotor Roo 96mm. 
Inner Radius of Outer Rotor Roi 80mm 
Outer Radius of Stator Rio 79.5mm 
Inner Radius of Stator Ri 42.5mm 
Outer Radius of Inner Rotor Roo 20mm 


The torque produced by the interaction between the it” PM 


pole and jt” coil can be written into the following expression: 
Tij = f(bij)di jij (1) 

$15 = cos 'ri8;/ |:|] lls; (2) 

d; j = r; x s;/|lri x s5|| (3) 


Where ¢;,; is defined as the included angle between the ith 
PM pole and jt” coil, f(¢;,;) is the torque constant function, 
r; is the position vector of the it” PM pole, and s; is the 
position vector of the j*” coil. T; j is the torque vector from 
the interaction of the it” PM pole and j*” coil. 
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Herein, the curve fit torque constant function f (Pij) has 
the following form: 


8 
f(¢i,3) = ao + ` (a; cos(ġi j * w) + bisin(pi j *w)) (A) 
i=1 
where ao, ai, b; and w are coefficients determined by the least 
square method | f (Pij) — F (Pij) | . The source function 


f(¢:,;) is obtained from the FE computed data as shown in 
Fig. 3. Parameters of the curve fitting equation are list on Table. 
II. The FE data comes from a pair of PM poles and one coil. 


e FETorque || 
= Fitting Curve} | 


Torque/Nmm 


0 5 10 1 20 235 30 35 40 45 
deg 


Fig. 3. Curve fit of the FE computed torque 


TABLE II. PARAMETERS OF THE CURVE FITTING EQUATION 
a0 0.5605 w | 0.04718 
al 0.6784 b1 | —0.6785 
a2 | —0.0841 | b2 | —0.7752 
a3 —0.4719 b3 —0.3074 
a4 | —0.3559 | b4 | 0.1259 


a5 | —0.09836 | b5 0.2108 
a6 0.01665 b6 0.1033 
aT 0.0133 b7 0.01937 


According to the symmetrical characteristics,when the coil 
is in different positions symmetrical about the PM axis, the 
torque between the coil and a pair of PM poles(from inner and 
outer rotors respectively) is of the same amplitude considering 
the situation of only one pair of PM poles and one coil. 
Spinning torque and tilting torque amplitudes are compared 
with the FE data. The spinning torque and tilting torque fits 
very well at the same included angle. The error range is less 
than 2%. Thus, through coordinate transformation, the torque 
between one coil and a pair of PM poles in the 3D space can 
be obtained. 


B. Torque Model of Single Coil in Rotor Frames 


Supposing (rx, "y, rz) is the direction vector of any axis in 
3D space. Then the rotation matrix R can be expressed as the 
Eqn. 5 when rotating a @ degree about the axis. With CEM, the 
analytical torque model can be obtained when the coil rotates 
about Z axis. FE results show that when the coil rotates in xy 
plane, the torque component in Z direction is much larger than 
in X and Y direction. Thus the torque model can be obtained 
from (0,0, Tz) or (0, Ty, 0) rotation. 


+ 


Fig. 4. Rotation of the Coil 


As shown in Fig. 4, spherical coordinate of the coil 
start point in xz plane is (r,0,0). The end point is 
(r, 0,0). Then, the torque vector can be obtained from the 
f (01) i (cos @1,0,sin@,) by rotating 92 about x axis. From 
Fig. 4, 01, 02 are the functions of 0 and ¢: 


6, = arccos(sin 6 cos p) (6) 
62 = arccos(sin 6 cos wy) (7) 


Herein,the form of torque rotation matrix is changed into : 


1 0 0 
Ro =| 0 coso sinə (8) 
0 —sinĝ cos 


The torque vector in (r, 0, ġ) can be calculated through Eqn. 
9 


T = |[Tz, Ty, Tz] = [0, F(61), Ro (9) 


As shown in Fig. 5 and 6, the maximum error between the 
modified analytical model and the FE data is less than 1%. 


Thus, the torque model of single coil in the rotor frames 
can be calculated through rotation and superposition. Through 
rotating for (i — 1)45°, the torque model between the i*” PM 
and coil can be obtained. It needs to be pointed out that 
the torque is opposite due to the alternate PM poles. The 
result shows that the torque produced by single coil has the 
mechanical circle of 90°. Fig. 7 and Fig. 8 show the results 
in the total rotating area. Where 6; represents the longitudinal 
angle, to illustrate the symmetry about the equatorial plane, 
90° has been deducted. Thus, 0 is from —20° to 20° and #1 
is from 0° to 360°. It can be concluded that the analytical 
model fits the FE results very well. The rotating torque is 
symmetrical about the equatorial plane, and is cyclical in the ¢@ 
direction. A total of four cycles exist in this direction. Because 
of the influences of the alternate PM poles, two peaks and two 
troughs exist in a cycle. 
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Fig. 5. Comparison of Analytical and FE Torque Results in y Direction 
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Fig. 6. Comparison of Analytical and FE Torque Results in z Direction 
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Fig. 7. Comparison of Analytical and FE Torque in z Direction Results 


(1—cos@)rery —rzsing (1—cos¢)rerz +rysing 
cos ¢ + (1 — cos bry 
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(1 — cos ¢)ryrz — rz sing (5) 
cos + (1 — cos ¢)r? 


Analytical 


Result FE Result 


T/Nmm 
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Fig. 8. Comparison of Analytical and FE Torque in Tilting Direction 


IV. FORMULATION OF MAGNETIC FORCE 


To benefit the improvement of the rotor support, magnetic 
levitation force produced by the interaction of PM poles and 
the coils can be utilized. Thus magnetic levitation force model 
is of obvious significance. In this section, it is discussed based 
on the FE simulation and CFM. The magnetic levitation force 
can be defined as the force exerted on the rotor in the radial 
direction by the coils. In other words, the levitation force 
produced by one coil is the radial projection of the counter- 
force from the energized coil. So our method is as following: 
first, the analytical force expression between one PM pole and 
one coil is obtained with FE data and CFM; second, according 
to the symmetrical characteristics, the model is extended to the 
whole 3D space. 


j; did © Si 8 


Frij ==gő (Yiz) j Is; | Is; | (10) 
J J 


Eqn. 10 can be departed into three parts: go (Wi,;) 2; is 
the amplitude of the force between the PM pole and the coil. 
go (Wi, j) ij is projection of the direction vector of the force on 
the radial direction. go is determined by Eqn. 11. The third 
part =n) is the direction vector of the magnetic levitation 


force. 


90 (Wig) = GH (11) 
O [Fiy 

S; = sj oR (12) 

cos(y) sin(y) 0 
R, = | —sin(y) cos(y) 0 (13) 

0 0 1 

al 
dij = RO Faia) Eya Faa) (14) 
[Fig 
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Fig. 9. Curve Fit of the FE Computed Torque 


Analytical Result 
FE Result 


bdeg * NOC” Adeg 


Fig. 10. Comparison of Analytical and FE Magnetic Leviation Force in x 
Direction Results 


Similar to the process of torque derivation,spherical coor- 
dinate of the coil start point in xz plane is (r, 0,0). The end 
point is (r, 0, œ). Then, the force vector can be obtained from 
the g (61) 7 (cos0,,0,sin@,) by rotating 92 about x axis. 61 
and 2 are defined in Eqn. 6 and 7. The final form of F,çi, J) 
can be written as: 

Fraj) = ij (Gxsin8 cos p + (Gycos02 — gz sin 02) sin 8 sin y 

+ (Gy sin 02 + gz cos 02) cos 0) (sin# cosy sinf sin y cos @) 
(15) 
The FE force data and the analytical curve of three components 
of Force in xz plane are shown in Fig. 9. The fit equation form 
is the same as Eqn. 4. The fourier fit type can satisfy the fitting 
precision very well. Fig. 10, 11 and 12 show the comparisons 
of analytical results with the FE results in 3D view. It can be 

concluded that the model is fit with the FE result very well. 


V. CONCLUSION 


A novel PM spherical actuator with three-dimensional pole 
array is proposed in this paper. The torque output and the 
magnetic leviation force of this spherical actuator are modeled 
based on the curve fit method with FE results. The analytical 


Analytical Result 
FE Result 


Fr_y/N 


Fig. 11. Comparison of Analytical and FE Magnetic Leviation Force in y 
Direction Results 
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Fig. 12. Comparison of Analytical and FE Magnetic Leviation Force in z 
Direction Results 


model mainly benefits the real-time control and improvement 
of the bearing of the spherical actuator. According to the 
geometric characteristics, the 2D fitted models of the torque 
and magnetic leviation force is successfully extended to the 
3D space. Results show good consistence of the analytical 
model and the FE results. The modeling errors are within 
the acceptance. Through appropriate correction, the modeling 
errors have been reduced to a lower lever. Thus, the analytical 
models can validate the FE results very well and can be used 
in the real-time control. 
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